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April 6, 2006. Note: some normalization errors were discovered in our computer code. In order 
to correctly interpret numerical results in this paper please see [lj. We investigate the effect of 
wiggly cosmic strings on the cosmic microwave background radiation anisotropy and matter power 
spectrum by modifying the string network model used by Albrecht et al. We employ the wiggly 
equation of state for strings and the one-scale model for the cosmological evolution of certain network 
characteristics. For the same choice of simulation parameters we compare the results with and 
without including wiggliness in the model and find that wiggliness together with the accompanying 
low string velocities lead to a significant peak in the microwave background anisotropy and to an 
enhancement in the matter power spectrum. For the cosmologies we have investigated (standard 
CDM, and, CDM plus cosmological constant), and within the limitations of our modeling of the 
string network, the anisotropy is in reasonable agreement with current observations but the COBE 
normalized amplitude of density perturbations is lower than what the data suggests. In the case of 
a cosmological constant and CDM model, a bias factor of about 2 is required. 



I. INTRODUCTION 

The origin of the large-scale structure of the Universe 
has been a focus of active research for the last few decades 
and is still one of the major unresolved problems in cos- 
mology. Measurements of the anisotropy in the Cosmic 
Microwave Background Radiation (CMBR) have inspired 
hope that different theories of structure formation can 
be tested obscrvationally. Quantum fluctuations in infla- 
tionary models and topological defects are two important 
sources of density inhomogeneities that would have dif- 
ferent imprints on the CMBR and which one could hope 
to distinguish. 

Cosmic strings are a special variety of topological de- 
fects that have been studied in the context of seeding 
large-scale structure. Rapid progress on both observa- 
tional and theoretical issues has resulted in more accu- 
rate estimates of the CMBR anisotropy due to strings 
[2-5,7,6]. The existence of publicly available computer 
codes [8] that can determine the anisotropy for given 
matter energy-momentum tensors has greatly facilitated 
the study of mechanisms that can source the CMBR 
anisotropy. 

The major hurdle in calculating the effect of strings on 
the CMBR is that there is no simple way to characterize 
the network of strings. Large scale computer simulations 
have provided insight into some of the properties of the 
network [9] within some choice of background cosmolo- 
gies. However, there are still unknown characteristics 
that could be important for predictions of the CMBR 
anisotropy and the large-scale power spectrum. As a re- 
sult, certain extrapolations have to be made to charac- 
terize the string network in the regime where details are 
not yet available. This "network modeling" is the crucial 
aspect of analyzing the observable signatures of cosmic 



strings. 

Over the last decade or so, a number of approaches 
have been taken to determine the detailed imprint of cos- 
mic strings on the CMBR. The first attempt in this direc- 
tion worked directly with the network of strings found in 
computer simulations [10]. Recently, this numerical anal- 
ysis has been refined and extended [11]. The drawback 
in any such numerical attempt is that the dynamic range 
of the simulations cannot yet extend to cover a cosmic 
expansion factor of at least 10 4 . It appears that some 
network modeling is essential as was first attempted in 
[12]. In recent work [5] the CMBR anisotropy was calcu- 
lated using the results of a lattice simulation of a string 
network in a flat spacetime background. Another ap- 
proach, first suggested in [13] and further developed in 
[3,4], approximates the string network by a collection of 
randomly oriented straight string segments, moving with 
random velocities. This model has the merit of being rel- 
atively simple and amenable to modifications that seem 
to be indicated by direct simulations. Hence, we have 
adopted this model and, for the first time, included small- 
scale structure on the string segments. Furthermore, we 
have modeled the parameters of the segments (length, 
velocity, wiggliness) by using the "one-scale model" of 
cosmic string networks [14,15] and some reasonable ex- 
pectations of how the strings would behave in an inflating 
background. 

The most important feature that we have taken into 
account in modeling the string network is the presence 
of small-scale structure on the strings. The scale of the 
wiggles is much smaller than the characteristic length of 
the string. In fact, a distant observer will not be able to 
resolve the details of the wiggly structure. Instead, s/he 
would see a smooth string, but the string would possess 
somewhat different characteristics. For a string with no 
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wiggles, the equation of state is fj, = T = /iq, where /i is 
the effective mass per unit length and T is the tension 
of the string. A wiggly string, however, has a different 
effective equation of state: [16,17], 

The wiggly string is heavier and slower than a Nambu- 
Goto string. 

The perturbation in the metric produced by a wig- 
gly string is similar to that of a smooth string. As was 
shown in [18], in both cases, the space in the vicinity of a 
straight segment of a string is conical with a deficit angle 
A = 8ttG[j,. The deficit angle is larger for wiggly strings, 
since they have a larger effective mass per unit length. 
However, a new important feature of the wiggly string 
metric is a non-zero Newtonian potential <fi oc (p, — T). 
(So the wiggly string behaves as a superposition of a mas- 
sive rod and a string with a conical deficit.) As a result, 
massive particles will experience a new attractive force 
F oc (p, — T)/v s , where v s is the string velocity [18]. At 
the same time, the propagation of photons in a plane 
perpendicular to the string is unaffected by the rod-like 
feature of the wiggly string metric and is only affected 
by the conical deficit. This fact - that the effect of wig- 
gliness on massive particles is qualitatively different from 
the effect on radiation - seems especially relevant when 
calculating CMBR anisotropy and the power spectrum. 
It may be hoped that this feature could help alleviate the 
current difficulties in reconciling the COBE normalized 
matter power spectrum with the observational data in 
the cosmic string model. 

It should be pointed out that the model we have 
adopted cannot be the final word since a number of fac- 
tors have still not been taken into account. In the net- 
work itself, the infinite strings will produce loops which 
will then decay by emitting gravitational radiation. The 
effect of the loops is only included insofar that they can 
be treated as segments of strings. In other words, large 
loops arc included in our analysis but small loops have 
been missed. Also, the decay of loops into gravitational 
radiation has been missed. The back-reaction of gravi- 
tational radiation on the string segments is completely 
neglected. We hope to rectify some of these omissions by 
further modeling of the network in subsequent work. 

In the following sections we use the model developed in 
Ref. [3] to calculate the CMBR anisotropy and the matter 
power spectrum due to cosmic strings in a few cosmolo- 
gies. The energy-momentum tensor that we use is that of 
wiggly strings, as we describe in the next section. In ad- 
dition, the characteristics of the string network (velocity 
and correlation length) are computed using the one-scale 
model described in Sec. III. There are some further as- 
sumptions that we have to make about the string network 
that we also describe in this section. Our results for the 
CMBR anisotropy and the power spectrum in density in- 
homogeneities for a few different cosmologies are given in 



Sec. IV. Here we also discuss the model dependence of 
our results. Our conclusions are summarized in Sec. V. 



II. WIGGLY STRINGS 

The world history of a string can be represented by a 
two-dimensional surface in spacetime, 

x" = x»(( a ), a = 0,1, 

called the string worldsheet, where £° and C 1 are the 
coordinates on the worldsheet. 

We consider strings that live in an expanding universe 
described by the metric, 

ds 2 =a 2 (T)(dr 2 -dx 2 ), 

where r is the conformal time. The string equa- 
tions of motion are invariant under reparametrization 
of the string worldsheet. A convenient choice of the 
parametrization (gauge) is 

C° = T , X' • X = 0, 

where the prime denotes a derivative with respect to C 1 . 
In this gauge the energy-momentum tensor of a string is 

T" v (y) = J d 2 ((ex»x v - e^x p x' v )5 {4) (y - x{Q), 

(1) 

where a — C} and 



e = 



(2) 



Following Carter [16] we can define the string tension 
T and the string mass-energy per unit length U by 

V=gT^(y) = J d 2 C,^(Uu»u v - Tv»v v )5 iA Hy - x{Q), 

(3) 

where and v p are such that 



UnW = -V..V 



= 1, u^ = 0, 



(u>*v p - v p u p )(u p v u - v p u u ) = <, 



(4) 
(5) 
(6) 



and r y ab is the worldsheet metric. One can check that 



(_ 7 )l/4< 



Vi(-7) 1/4 



(7) 



satisfy the conditions (4)-(6). Substituting (7) into (1) 
and comparing with (3) one obtains 
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U = T = n, 



(8) 



which is the equation of state for an ordinary ( "intrinsi- 
cally isotropic" [16]) Nambu-Goto string [19]. 

Lattice simulations of string formation and evolution 
suggest that strings are not straight but have wiggles 
[9]. A distant observer, however, would not be able to 
resolve the small-scale structure. Insead, he would see a 
smooth string with the effective mass per unit length U 
and the tension T. The equation of state for a wiggly 
string averaged over the small-scale structure has been 
shown [16] [17] to be 



or 



UT = /i 2 



U = cefl, T = [l/oi, 



(9) 



(10) 



where a — a(a,r) is, in general, some function of time 
and the coordinate along the length of the smoothed 
string. The energy-momentum tensor of a wiggly string 
viewed by an observer who cannot resolve the wiggly 
structure is then obtained by substituting (10) into (3): 

f^(y) = -j= J ^=.{Uu»u v - TW)<^(y - x(Q) 



V=9 J 



Next, we use this expression to calculated the stress- 
energy of a string network. 



III. THE STRING NETWORK 

A. Parameters of the network 

The string network at any time can be characterized 
by a single length scale, the correlation length L, defined 

by 



L 2: 



(11) 



where p is the energy density in the string network. It 
will be convenient to work with the comoving correlation 
length I = L/a. 

The expansion stretches the strings, thus increasing the 
energy density. At the same time, long strings reconnnect 
and chop off loops which later decay. The evolution of / 
including only these two competing processes is [14,15,20] 



dl a o 1 
— = -Iv + -cv 
dr a 2 



(12) 



(13) 



where v is the rms string velocity, c is the loop chopping 
efficiency and k is the effective curvature of the strings. 
The values of c and k in radiation and matter eras are 
suggested in [20] . We use the same scheme as in Ref. [4] to 
interpolate between these values through the radiation- 
matter transition: 

c r + gac m 



c(r) 
k(r) = 



l+ga 
+ gak m 



l+ga 

where we take cy = 0.23, c m = 0.18, k r = 0.17, 
k m = 0.49, g = 300 and a(r) is normalized so that a = 1 
today* . 

B. Model of the network 

Here we are closely following the model described in 
[3,4]. 

The basic picture is that the string network is repre- 
sented by a collection of uncorrelated, straight string seg- 
ments moving with random, uncorrelated velocities. All 
the segments are assumed to be produced at some early 
epoch. At every subsequent epoch, a certain fraction of 
the number of segments decay in a way that maintains 
network scaling. This picture of the network is depicted 
in Fig. 1. 

The comoving length, I, of each segment at any time 
is taken to be equal to the correlation length of the net- 
work defined below eq. (11). The positions of the seg- 
ments are drawn from a uniform distribution in space 
and their orientations are chosen from a uniform distri- 
bution on a two sphere. The segment speeds are fixed to 
be given by the solution of eq. (13) while the direction 
of the velocity is taken to be uniformly distributed in the 
plane perpendicular to the string orientation^. In princi- 
ple, this constraint on the velocity does not remain valid 
when the strings are wiggly since the wiggles can impart 
a longitudinal velocity to the segments. However, as ex- 
plained below, the longitudinal velocities are expected to 
be much smaller than the transverse velocities and hence 
will be neglected. 

The decay of segments of the string network is accom- 
plished by "turning off" the energy-momentum of a frac- 
tion of the existing segments at every epoch. Each seg- 
ment is assigned a certain decay time, r m , where the 



*This choice of the value of g along with our normalization 
for a leads to the same time-dependence of v and I as in [4]. 
Normalizing a so that a = 1 at equality would require a much 
smaller value, e.g. g ~ 0.1, as reported in [4] 

*We have also performed a few simulations where we drew 
the velocities from a Gaussian distribution as in Ref. [3], but 
these did not lead to significantly different results. 
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index m labels the individual segments. So the Fourier 
transform of the total stress-energy of the network is the 
sum over the stress-energies of all segements: 



e^(k, t) = £ (k, t)T oS (t, r m ) , (14) 

m— 1 

where -/Vo is the initial number of segments, and 
T off (r, r m ) is a smooth function that turns off the m th 
string segment by time r m * . The functional form is taken 
to be [3] 



T off (r,r m ) 




3x) 



. t < Lr m 
. Lr m < t < 

■T > T m 



where 



HLTjn/r) _ 

ML) 



(15) 



(16) 



and L < 1 is a parameter that controls how fast the 
segments decay. 

The total energy of the string network in a volume V 
at any time is 



N/j,L = Vp 



L 2 



(17) 



where N = N(t) is the total number of string segments 
at that time, V = Vba 3 , a = 1 at the present epoch and 
Vo is a constant simulation volume. From (17) it follows 
that 



L 3 Z 3 • 



(18) 



The comoving length I is approximately proportional 
to the conformal time r and implies that the number 
of strings N(t) within the simulation volume V) falls as 
t -3 . To calculate the CMB anisotropy we need to evolve 
the string network over at least four orders of magnitude 
in cosmic expansion. Hence we would have to start with 
N > 10 12 string segments in order to have one segement 
left at the present time. This is the main problem in 
directly dealing with the expression for the energy mo- 
mentum tensor in eq. (14). 

A way around this difficulty was suggested in Ref. [3] . 
The suggestion is to consolidate all string segments that 
decay at the same epoch. Since the number of segments 
that decay at the (discretized) conformal time tj is 



*In addition, Albrecht et al. introduce a function T on which 
turns on the segments at a very early time. This function is 
not a feature of the model but only introduced to speed up 
the code. We have not included T on in our simulations. 



string segments 

FIG. 1. A schematic picture of the string network model. 
All string segments (depicted by hlled circles) are born at 
an early epoch and then decay at various later times. The 
segments are labeled by the index m and the decay times are 
shown as numbers along the r axis. In certain cosmologies, it 
is possible that some string segments never decay. 



consolidated string segments 

FIG. 2. The modified model of the string network. All 
strings that decay at the same (discretized) time in Fig. 1 are 
consolidated into one string segment and assigned a weight 
that is the square root of the number of segments that the 
consolidated segment represents. This works for all segments 
that decay by the end of the simulation but will miss those 
segments that do not decay. In the present scheme, the seg- 
ments that will not have decayed by the end of the simulation 
are consolidated into one string. The contribution of this sur- 
viving segment is the remainder term in the expression for 
the energy momentum tensor in eq. (20). 
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N d (n) = V[n(T^ 1 )-n(n)} 



(19) 



where n(r) is the number density of strings at time r, 
the "consolidated string" decaying at n is taken to have 
weight: y/Ndin). (The assumption is that string seg- 
ments that decay at the same time act randomly, leading 
to the square root in the weight.) The number of consol- 
idated strings is of the order of a few hundred and can 
be dealt with computationally^ . This modified picture of 
the string network is shown in Fig. 2. 

Now we can choose r„ to be equally spaced on a log- 
arithmic scale between r m i n and r max and write the 
energy-momentum tensor as 

K 
i=l 

(20) 

where K is the number of consolidated segments and 
is a remainder that we have included and that we now 
explain. 

The sum in eq. (20) misses any string segments that 
have not decayed by T max . The remainder is sup- 
posed to represent the contribution of these segments 
of strings. As in the case of the decaying strings, we 
will consolidate these surviving segment into one segment 



with the weight \JVn(T max ). Therefore 



V = VVn( Tmax ) 0$(fc,T) 



(21) 



where, the superscript on 6 means that it is the energy- 
momentum tensor for one string segment. 

For certain cosmologies - those without a cosmological 
constant - the remainder term can be made arbitrarily 
small by choosing a large enough T max . This is because 
n(r) oc Z(t)~ 3 , where 1(t) <~ r is the comoving correlation 
length, and hence n(r) — > as r — > oo. However, if the 
cosmological constant is non-zero, the universe enters an 
inflationary epoch and the range of the conformal time 
is finite, r € [0, Too], an d given by: 



dr 



dt 

W) 



In this model, the number density of string segments 
is assumed to only depend on the correlation length as 
n{r) oc Z(t)~ 3 . However, the decay of string segments as 
described by the function prevents this relationship 
from being exact. Instead we have: 



n(r) = 



C{r) 
l{rf 



5 In addition to consolidating the string segements decaying 
at a given time into a single "consolidated string" we have 
also tried consolidating them into two, three and four strings. 
This did not lead to any difference in the final results. 



The function C(t) is determined by requiring that the to- 
tal number of strings at any time is still given by V/1(t) 3 . 
Therefore: 



1 K 

-Ts = J^nfa-i) -n(T 4 )]T off (T,T 4 ) + n(T maa; ), 



and l(r) is determined from the one-scale model (eq. 
(12)). (Note that C(r) is contained in n{r).) It is reas- 
suring to see that the function C (t) that we obtain from 
our code is nearly constant and of order unity throughout 
the simulation. 

The Fourier transform of the energy- momentum tensor 
of an individual string segment is 



e„„(fe,T) = J 



d 6 xe lkx Q^ 

1/2 



1 



H / dae ik - x UaX^X u X'^X W ). 

-1/2 «* 



(22) 



where, X^{a,r) are the coordinates of the segment: 



X 1 



t, X = x + crX' + vtX, 



(23) 



where xo is the random location of the center of mass, 
X' and X are randomly oriented unit vectors satisfying 

X' ■ X = and v is the velocity of the string. 

The random location vectors, Xq, appear only in the 
dot product Xo ■ k. Instead of generating x at random, 
Albrecht et al. generate xq ■ k for each segment as a 
random number from [0, 27r]. We have also followed this 
scheme**. 

As mentioned earlier, we will be ignoring the longitu- 
dinal velocities of the segments. This is the constraint 
X' ■ X = that we have imposed. The justification is 
that there are on the order of 10 4 left- and right- mov- 
ing wiggles on a segment of string, each traveling with a 
velocity of about the speed of light. The average momen- 
tum of this "gas" of wiggles will be zero but there can 
be fluctuations which will lead to a velocity ~ 1/VlO 4 . 
Such longitudinal velocities are insignificant compared to 
the transverse velocity expected to be ~ 0.1. 

The differential equations which describe the metric 
perurbations produced by 6^(fc, r) do not depened on 



**In ongoing work, one improvement of the model that we 
are investigating is a scheme in which we choose k m i n xo ■ k at 
random in the interval [0, 2n] and then multiply by k/k m i n 
where k min is the smallest wave-vector considered in the 
simulation. 
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the direction of k. Therefore, without any loss of gener- 
ality we will assume k — k^k in cq. (22). 

Substituting (23) into (22), integrating over a and tak- 
ing the real part gives 



e 00 



lia sin(/cX^/2) 



Vi^2 

v 2 XiXi 



kX'j2 
(1-t, 2 ) 



cos(k ■ x + kX 3 vr) , (24) 
6oo. (25) 



2 X ' X J 



Ooi can be found from the covariant conservation equa- 
tions VS^ = 0. Here we have assumed that the entire 
string energy is in long strings. The conservation is vio- 
lated if one includes the energy in loops and the gravita- 
tional radiation. A detailed discussion of this issue can 
be found in [5]. 

C. Evolution of parameters 




The evolution of the velocity, v, with time is found 
by solving eqns. (12) and (13) with a range of possible 
intitial conditions. In Fig. 3 we show the evolution of 
the string velocity for the initial conditions: l{T m in) — 
0.13r m i„, vo(T m i n ) — 0.65 for two different cosmological 
models. Similarly we show the behaviour of 1(t)/t as a 
function of r in Fig. 4. 

The values of the "wiggliness" parameter a have been 
estimated [9] in the radiation and matter eras to be a r ~ 
1.9 and a m ~ 1.5 respectively. In the case of a non-zero 
cosmological constant, the exponential expansion of the 
Universe will stretch and smooth out the wiggles, so that 
a — > 1 in the A— dominated epoch. A function that fits 
the expected evolution of a is (Fig. 5) 



a(r) = 1 + 



(a r — l)a 



(26) 



We shall assume this behaviour of a in the next section. 

Next, we use the expression for the string stress-energy 
to calculate the radiation and matter perturbation spec- 
tra at present time. 



FIG. 3. The velocity of strings v as a function of the confor- 
mal time for flbaryons = -05, Ocdm = 0.95, Qa = (dotted 
line) and fltaryons = -05, 0.cdm = 0.25, Qa = 0.7 (solid line). 



IV. RESULTS 



We calculate the CMB anisotropy using the line of 
sight integration approach [8] implemented in the pub- 
licly available code CMB FAST. Given 9^^ for each k and 
r, CMBFAST integrates the Einstein equations simulta- 
neously with the Boltzmann equations for all radiation 
and matter particles present in the model. Two factors 
that have not been included in our analysis are: (1) com- 
pensation of string inhomogcncities by perturbations in 
the matter fluids, and, (2) small loops of string. The first 
factor was considered in Ref. [3] which we are closely fol- 
lowing. However, it was found there that the effects of 



1 10 100 1000 io 4 



FIG. 4. The length of string segments, l(r), divided by r 
as a function of r for Qbaryons = -05, Sic dm = 0.95, J1a = 
(dotted line) and flbaryons = -05, 0.cdm = 0.25, Qa = 0.7 
(solid line). 
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FIG. 5. The wiggliness a as a function of the conformal 
time for tlbaryons = -05, 0.cdm = 0.95, = (dotted line) 
and Q,bar V ons = -05, Hcdm = 0.25, £Ia = 0.7 (solid line). 



compensation do not significantly alter the predictions of 
anisotropy and power spectrum. It is possible that the 
second factor is important. At the moment, we do not 
have a good model for including loops in our analysis as 
there is not much information available in the literature 
on which to build such a model. The inclusion of loops 
in the model is a problem that we postpone for the fu- 
ture, though the slow decay of string segments described 
below may mimic the presence of loops to some extent. 

The CMBR anisotropy as seen by an observer can be 
described by A(x, n, r Q ) = |T(x, h, t q ) — T\/T, where x is 
the position of the observer, n is the line of sight direc- 
tion, t is the conformal time today and T is the average 
temperature of CMBR. At x = one can decompose A 
into spherical harmonics 



A(n) = ^2ai m Y hn (n). 



(27) 



Irn 



The angular power spectrum Ci is denned by 



Q = 



l 



21 + 1 



X] ( a lm a lm), 



(28) 



where () denotes an ensemble average. A self-consistent 
treatment of the line of sight method and the complete 
set of the equations is given in [21]. The output of CMB- 
FAST is the COBE normalized angular power spectrum 
of temperature anisotropy and the power spectra of each 
of the matter species that is present. 

In its current version CMBFAST calculates perturba- 
tions from scalar and tensor sources. Since strings may 
generate a significant vector component, we had to add 



1 1 1 
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i i 1 
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FIG. 6. The total angular power spectrum with (solid 
line) and without (dashed line) including the wiggliness when 
tibaryons = -05, Q,cdm = 0.95, Q.A = and vo = 0.65. The 
compiled observational data is plotted for comparison. 



j ' ' I 1 



FIG. 7. The matter power spectrum for the same choice of 
parameters as in Fig. 6. The data extracted from surveys of 
galaxies and clusters of galaxies is plotted for comparison. 
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the vector part to the code^. 

Our results for the anisotropy and power spectrum are 
averages over M different realizations of a random net- 
work of strings. In the limit of large M, each quantity 
that we calculate will have an associated probability dis- 
tribution. For example, for each I the distribution of G 
will approach some probability function with mean Gi 
and standard deviation 07. For the averaged result, Ci, 
the standard deviation is aij\[M. Performing a large 
number of experiments on a computer allows us to find 
Ci and 07 quite accurately. However, an observer (eg. 
COBE) trying to determine Ci would only be able to av- 
erage over (21 + 1) causally disconnected patches on the 
sky. Therefore, the accuracy in determining (7; observa- 
tionally is + 1. We take this limitation (known 

as "cosmic variance") into account when plotting the 
1 (j— error bars on our graphs. 

Our results were found by averaging over M = 300 
string network realizations, however, 100 runs would be 
enough to reproduce most of our results. The simulation 
with 100 runs took about 10 hours to run on a single 
processor IBM 590 workstation and about 5 hours on 
10 processors of the J90 cluster of the National Energy 
Research Scientific Computing Center. 

We work only with a flat universe with and with- 
out a cosmological constant. We consider two cases: 

(Qbaryons = -05, VLqDM = 0.95, f^A = 0) and (Qbaryons = 

.05, flcDM = 0.25, f^A = 0.7). The first case is what un- 
til recently was known as the standard CDM model. The 
second is motivated by the recent supernovae data that 
suggest Oa = 0.7 and a flat universe. The value of the 
Hubble constant in both cases was taken to be Hq = 50 
km sec^ 1 Mpc^ 1 . 

In Fig. 6 we show the effect that wiggliness has on 
C/'s for the case when tt\ = along with the compiled 
experimental data [22]. The matter power spectrum for 
the same set of paramters is shown in Fig. 7 together 
with the data from surveys of galaxies and clusters of 
galaxies [23]. Adding wiggliness results in a higher peak 
at I i=a 400. However, it does not improve the shape 
or the magnitude of the matter power spectrum which 
appears to be in disagreement with data. The COBE 
normalization allows us to determine the string mass per 
unit length. We obtain that G^o = 1-1 x 10~ 6 for wiggly 
strings and G[io = 1-5 x 10~ 6 for smooth strings. 

Allowing for a non-zero cosmological constant im- 
proves the agreement with the data. In Fig. 8 and Fig. 
9 we plot the results of the simulation with J7a = 0.7. 
The peak in the CMBR anisotropy is now higher com- 
pared to the il\ = case. The power spectrum in Fig. 9 
includes the factor of r^a tter necessary for correct com- 



"Our code can be downloaded from the website: 
http://theory4.phys. cwru . edu / ~levon . 
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FIG. 8. The total angular power spectrum for wiggly (solid 
line) and smooth (dashed line) strings when tlbaryons = -05, 
SIcdm = 0.25, Q.a = 0.7 and v = 0.65. 

parison with data [23]. Although the magnitude is sig- 
nificantly lower than the data, the shape of the matter 
power spectrum is improved. The wiggliness results in a 
higher peak in C;'s and a slight increase in the magni- 
tude of matter power spectrum at smaller length scales. 
The string mass per unit length obtained from COBE 
normalization is practically the same as in f^A = case. 

There is evidence from simulations [9] that the large 
scale string velocities are much smaller than those plotted 
in Fig. 3. Also, wiggly strings must be heavier and slower 
than smooth strings because of the different equation of 
state. We have modified the parameters c r , c m , k r and 
k m of the one-scale model in such a way that the rms 
velocity is 0.12 in the radiation era and 0.1 in the matter 
era. In Fig. 10 and Fig. 11 we plot the results of using 
smaller velocities. We see a significant increase in the 
magnitude of the matter power spectrum. Including the 
wiggliness results in a slight shift of power to smaller 
scales and an overall increase in magnitude. The scale- 
dependent bias, required to fit the data varies between 
b w 1.6 on smaller length scales and b s=s 2.4 on larger 
scales. COBE normalization gives G^o ~ 2.3 x 10~ 6 for 
smooth strings and G/jq ~ 1-9 x 10~ 6 for wiggly strings. 

In addition, we have tested the dependence on two 
other parameters of the model that may have physical 
significance. These are the string decay rate L and the 
correlation length of the network 1(t) . The parameter L 
appears in the function T°tf, eq.(15), and controls how 
early and how fast the string segments decay. Smaller 
values of L correspond to an earlier and slower decay. We 
have used L = 0.5 for most of our calculations (Figures 6- 
11,14,15). In Fig. 12 and Fig. 13 we compare the results 
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FIG. 9. The matter power spectrum for the same choice of 
parameters as in Fig. 8. 



of using L = 0.1, 0.5 and 0.92. Lower values of L tend to 
raise the matter power spectrum. When string segments 
start to decay earlier the effective string density at all 
times (and especially at later times) will decrease and a 
larger COBE normalization factor is required to fit the 
angular power spectrum on large scales. This, of course, 
encreases the value of G/j,q. We obtain G/io ~ 4 x 10~ 6 
when using L = 0.1. 

The evolution of the correlation length 1(t) (cq. (12)) 
is completely determined by the parameters of the one- 
scale model c r , c m , k r and k m . To obtain the results 
shown in Figs. 6-13 we have chosen these parameters 
so that the evolution of 1(t) is not altered and is the 
same as in Fig. 4. In Fig. 14 and Fig. 15 we plot 
the results of the simulation in which the cosmological 
parameters and the string velocity were taken to be those 
in Fig. 10 but the parameters of the one-scale model 
were set to give 1(t)/t equal to 0.24 in the radiation- 
and 0.3 in the matter-dominated eras. This choice of 
parameters is largely arbitrary, however, it shows that 
there is a freedom in the model that may be used to 
change the shape of the matter power spectrum. 



V. CONCLUSIONS 

We have seen that including the effects of small-scale 
structure in the string stress-energy improves the agree- 
ment with the observational data. Using ft a — 0.7 and 
reducing the average string velocities gives us a scale- 
dependent bias factor between 1.6 and 2.4. Allowing 
for a slower decay rate of strings and larger correlation 
length of the string network make the bias factor nearly 
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FIG. 10. The total angular power spectrum for wig- 
gly (solid line) and smooth (dashed line) strings when 
flbaryons = -05, flcDM = 0.25 and f^A = 0.7 and using small 
values for the string velocities: v = 0.12 in the radiation era 
and 0.1 in the matter era. 




FIG. 11. The matter power spectrum for the same choice 
of parameters as in Fig. 10. 
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FIG. 12. The total angular power spectrum for the net- 
work of wiggly strings with L = 0.1 (dotted line), L — 0.5 
(long dash line) and L = 0.92 (short dash line). All other 
parameters are the same as in Fig. 10. 

scale- invariant with an average value of 1.9. The scale- 
invariant bias factor of 2 was reported in [4] which was 
obtained by introducing an additional time-dependence 
of n designed to reproduce the shape of the matter power 
spectrum. The present work has the advantage of not al- 
tering the physics of the model. 

The one-scale model, eqns. (12) and (13), is a very 
coarse approximation of the evolution of wiggly strings. 
The small-scale structure is included only through the 
value of the effective curvature k, the rest of the equa- 
tions being exactly the same as for the smooth strings. 
We expect that a more precise analytical model would 
improve the agreement with the data. Also, work on im- 
proving the model of the string network is in progress 
as reported in [5] and this should help to eliminate some 
of the assumptions we have had to make in the present 
analysis. In addition we are presently investigating ways 
to relax some of the assumptions of the model described 
in the paper that will help shed light on the robustness 
of the predictions. 
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FIG. 13. The matter power spectrum for the same choice 
of parameters as in Fig. 12. 
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FIG. 14. The total angular power spectrum for the network 
of wiggly strings with L = 0.1. The solid line corresponds to 
the modified 1(t). The dotted line is the same as in Fig. 13. 
Other parameters are the same as in Fig. 10. 
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FIG. 15. The matter power spectrum for the same choice 
of parameters as in Fig. 14. 
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